Disorder-Driven Magnetic Field Dependence of the Internal Field Distribution in the 

Bragg Glass Phase of Type-II Superconductors 
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We use the replica variational method to study the effects of weak point disorder on the variance 
of the internal field distribution measured in NMR and muon-spin rotation experiments in type-ll 
superconductors. We show that for a simple model there is significant magnetic field dependence 
which is extrinsic and disorder-driven, and does not have a microscopic (non s— wave pairing) origin. 
Results are presented where we examine the dependence of the magnetic field variance upon the 
strength of the applied external field. 
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The study of high temperature superconductors (high- 
Tc's) remains a subject of intense investigation, both in 
the search of an explanation of its microscopic origin and 
in the application of phenomenological theories to investi- 
gate the effects of thermal fluctuations, random disorder, 
and transport. It is now generally accepted that high-Tc 
superconductors are unconventional d— wave rather than 
conventional s— wave with an energy gap that vanishes 
along certain directions in momentum space, resulting 
in nodes in the superconducting gap. Measurement of 
the temperature dependence of the magnetic penetra- 
tion depth, A(r), is one way to probe the non s— wave 
nature of the low energy nodal excitations of the super- 
conducting state [1]. A detailed understanding of the 
d— wave nature of the high-Tc's must also include a quan- 
titative description of the dependence of the two fun- 
damental length scales characterizing a superconductor, 
the penetration depth A and the coherence length ^o, on 
the strength of an applied magnetic field Bq. An early 
theoretical investigation of the weak-field response of a 
da;2_y2 superconductor predicted a direction-dependent 
nonlinear Meissner effect, associated with the quasiclas- 
sical shift of the excitation spectrum due to the super- 
flow created by the screening currents [2] . More recently, 
Amin et al. studied nonlinear and nonlocal effects due 
to the field induced excitations at the gap nodes to pre- 
dict the temperature and field dependence of an effective 
penetration depth [3]. 

While understanding the microscopic physics at play in 
the high-Tc's is a formidable challenge, both experimental 
and theoretical investigations have shown that the mag- 
netic field-temperature {B — T) phase diagram of these 
materials is also quite a bit more complicated than for 
conventional superconductors [4]. It is well known that 
the flux lines in a clean s— wave type-II superconductor 
form an Abrikosov triangular lattice for fields B larger 
than the lower critical field, B^ ■ However, for high-Tc's, 
this mean field phase diagram is considerably affected by 
the combined effects of thermal fluctuations and various 
types of disorder, such as oxygen vacancies, columnar 
defects, or twin boundaries (see Fig. la) [4]. 
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Fig. 1: (a) Schematic phase diagram of high temperature 
superconductors with point disorder, (b) The generic behav- 
ior of the relative displacements correlation function. 

Thermal fluctuations lead to melting of the vortex lat- 
tice and appearance of a vortex liquid (VL). In addi- 
tion, extensive experimental, theoretical and numerical 
work strongly suggest that the Abrikosov vortex lattice 
in type-II superconductors is destroyed by the presence 
of disorder and gives place to various superconducting 
glass- like mixed states [4]. If randomness is strong, the 
underlying translational order of the vortex lattice is 
completely destroyed giving a vortex glass phase (VG). 
However, for weak enough randomness, 'some order' is 
expected to survive. Indeed, it has been suggested that 
there may exist a stable, dislocation-free Bragg glass 
(BG) phase at low magnetic fields and temperatures in 
the presence of weak point impurities where the struc- 
ture factor 5(k) has power law singularities at the Bragg 
reciprocal vectors k = G [5,6]. For an applied field larger 
than a critical value B* (of order 2 to 20 Teslas, depend- 
ing on the specific sample considered) [7], the system 
undergoes a transition to a vortex glass which consists of 
completely randomly frozen vortex lines with finite width 
peaks in the structure factor at k « G (see Fig. la) 
[5,6,8]. 

The vortex state of high-Tc's leads to a spatially vary- 
ing magnetic field B(r) inside the superconductor which 
may be investigating using various techniques. The muon 
spin rotation (/iSR) technique is an almost ideal method 
for studying the magnetic field distribution in the mixed 
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state, and magnetic field and temperature dependence 
of A(i?o) and ^0(^0 ) [9]- In a /iSR experiment the in- 
homogeneous magnetic field is sensitively monitored by 
implanting muons into the sample, one at a time, where 
the muon's spin precess with a Larmor frequency which 
is directly proportional to the local magnetic field B(r). 
The muon-spin precession signal, which is obtained by 
detecting the muons positron-decay pattern, reflects the 
distribution of the Larmor precession frequencies from 
all the muon stopping sites, hence forming an asymmet- 
ric /iSR lineshape [9]. The width of this lineshape and, 
correspondingly, the width of the internal field distribu- 
tion, is roughly proportional to A parametrization 
of the experimental lineshape based on a theoretical line- 
shape produced by an effective London model that in- 
cludes both the cut-off effects arising from the finite-size 
of the vortex core and a phenomenological description 
of the effects of broadening caused by static and ran- 
dom distortion of the vortex lattice and/or nuclear dipo- 
lar magnetic fields allows for a determination of effective 
X^'^iBo) and (Bo) [9-11]. 

/uSR experiments on YBa2Cu306-i-5(YBCO) [9,11] 
have reported a magnetic field Bq dependence of the /iSR 
lineshape. This has been so far attributed to a micro- 
scopic and intrinsic magnetic field dependence of the Lon- 
don penetration length A and coherence length ^0 that 
could possibly be due to the underlying d— wave symme- 
try of the superconducting pairing state of YBCO [3]. 
This brings up the following question: Could the mag- 
netic field dependence of the /xSR lineshape and the field 
dependence ascribed to the extracted effective A°^(i3o) 
and Q^{Bo) be partially extrinsic and unrelated to the 
underlying microscopic rf— wave physics, but rather due 
to the fact that the lineshape ought to change from 
(somewhat) asymmetric in the Bragg glass state to more 
or less a Gaussian shape in the vortex glass phase? In- 
deed, such field-driven evolution of the asymmetry (re- 
ferred to as skewness) of the /xSR lineshape has been 
reported for BSCCO [12,13] and YBCOg.eo [14] high-T^ 
materials. In this paper, we study the applied magnetic 
field So dependence of the variance of the field distri- 
bution, and consequently A°^(_Bo), caused by weak point 
disorder. To this end, we use replica variational solutions 
of the displacement correlations of the flux lines in the 
presence of weak point disorder, and calculate the mag- 
netic field variance, which is related to the structure fac- 
tor S'(k) of the vortex system. We note, however, that 
the vortex-vortex displacement correlation functions of 
the weakly disordered vortex lattice is anisotropic and 
may have nonuniversal exponents [15]. Consequently, we 
use an isotropic model to simplify calculations and ex- 
pose the basic underlying physics at stake. Here, we only 
study the effect of disorder; the role of thermal fluctua- 
tion effects will be discussed elsewhere. 

The influence of the randomness on the translational 
correlation function has been discussed by many authors 
[4]. It was originally discussed by Larkin [16] using a 
model in which weak random forces act independently 



on each vortex. This model predicts an exponential de- 
cay of the translation correlations on length scales larger 
than a disorder dependent Larkin length, Rl- However, 
Bouchaud et al. pointed out that at large scales the lat- 
tice starts behaving collectively as an elastic manifold in 
a random potential with many metastable states and that 
the exponential decay does not hold beyond the Larkin 
length Rl [17]. They used a variational replica field the- 
ory to study the pinning of vortex lattice by impurities 
they found a power law roughening of the lattice with 
stretched exponential decay of the translational correla- 
tion fimction for length scale beyond Rl [17]. Later, 
Giamarchi and Le Doussal [5] showed that while dis- 
order produce algebraic growth of displacement corre- 
lations at short scales, periodicity takes over at large 
scales and results in a growth of displacements that is 
at most logarithmic in x. The results for the displace- 
ment correlation function are summarized in Fig. 1 (b). 
For length scales smaller than Rl the model is equivalent 
to Larkin's model with correlations cx x (in three dimen- 
sion) and which corresponds to the replica symmetric 
part of the variational solution. For Rl < x < Rc, the 
system is in the random manifold (RM) regime [17]. In 
this regime, replica symmetry is broken due the various 
metastable vortex configuration states and relative dis- 
placements correlate as x"^", where 1/ = 1/6 in variational 
approximation. For x > Rc, the periodicity of the lattice 
becomes important and one enters the asymptotic Bragg 
glass logarithmic regime [5,6]. 

In the London limit, the field of point vortices sitting at 
arbitrary positions r„j = {x„i, jjrm Zm) for A >> rf, where 
d is the distance between layers, varies slowly between 
layers, and is given by [18-20]: 

--3b(k)^exp[»k.(r-r„)] 

^ ^ m 

^^^>- fci(l + A2F) ' 

where k = (kj_, fc^;). 

Brandt [20] had used Eq. (1) to study the magnetic 

field variance in layered superconductors. Brandt showed 
that perturbation of a lattice of rather stiff flux lines in- 
crease the field variance [21]. However, the fluctuations of 
vortex line segments or vortex pancakc^s of highly fiexible 
flux lines may decrease the /iSR linewidth [20,22]. 

The general expression for the magnetic field variance, 
= ([B(r) - [B(r)]]2), where [...], < ... >, and ^ de- 
note space, thermal, and disorder average, respectively, 
can be written as 




(2) 

where V is the volume of the sample. The lattice sum 
is the structure factor of the vortex-point arrangement. 
One can rewrite Eq. (2) as 
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5^ = 8^ f^l b(k) r S(k), 



d J (27r)3 



(3) 



where po = -So/0o is the density of vortices and the struc- 
ture factor ^(k) is given by 

5(k) = (^) E/ '^'^ . (4) 

G 

Here G is reciprocal lattice vector and Cg{x) is the trans- 
lational correlation function, which for simple isotropic 
case has the following form [5,6] 



CG(x)=exp|^([«(a;)-«(0)]2>| , 



(5) 



where u{x) is the displacement of the flux line from its 
equilibrium position. 

The above equation for (5^ has in general both perpen- 
dicular, (5_L, and z, S^, components. Before analyzing the 
efii'ect of impurities on the magnetic field variance, we first 
discuss a simple limiting case to estimate the order of dif- 
ferent components. If the point vortices in each layer are 
assumed to be randomly positioned, then one simply has 
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(6) 



The first term gives the fluctuations of B^, and the sec- 
ond term the fluctuations of B_l. A cutoff has been 
used due to the factor for which the Joscphson 

length A J — dXc/X was inserted by replacing by 
l/{kj_ + [20]. Integrating over k leads to 



Bod(j}o 
87rA3 



cos-HAj/A) 
V(VA7)^ - 1 



(7) 



For Aj << A the second term which comes from the ran- 
dom perpendicular fluctuations of B(r) is small. Conse- 
quently, and also for sake of simplicity, the perpendicular 
contribution shall be disregarded in the rest of the paper. 

Returning to the general expression Eq. (2), using Eqs. 
(3), (4) and (5), and considering only the z component, 
and integrating first over all (six) angles and k we find 
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xe ''/^CG{x)sm{Gx)dx. (8) 



The above equation can then be used to study the vari- 
ance of the magnetic field. The simple case of Cg{x) = 1 
corresponds to the zero disorder case, in which inte- 
grating over X leads to the London model with 5'^ = 
_Bq X]g#o 1/(1 + A^G^)^. For finite disorder, recalls that 
the relative displacement correlation function has differ- 
ent behavior in the Larkin, RM, and logarithmic regimes 
(Fig. lb). Therefore, we have to carry the integration of 
Eq. (8) for these three x— dependent regime of Cg{x). 
Using expression (9) of Ref. [6] for Rc, 



Rc 



2„4 3/2 1/2 
^"0'-66 '-44 

7T^plU^2TTd^0 ' 



(9) 



where oq is the lattice spacing, is superconducting 
coherence length. Up is a typical pinning energy per 
unit length along z, and C44 and cee are elastic con- 
stants. Considering C44 = c£o/(7^aQ) and cge = £o/(4ao) 
with the vortex line tension Sq = ((/>o/47rA)2, and c = 
ln(7d/^o), 7 is the anisotropy ratio, one gets 



Rc — 



2cV2 



(10) 



The Larkin length is also given by Rl ~ {^o/ctoYRc- 
Taking parameters for YBCO 7 = 5, A = 1200A, 

d = IQA, ^0 = 10^, and Up/eo = 1/40, we see that 
Rl « 5l and Rc ~ lO^A (for Bq ~ IT). Since 
Rl « do, Rc » do, and Rc » A, the effects 
of Larkin and logarithmic regimes are very small, the 
main contribution to the integral Eq. (8) comes from 
the RM regime as we have confirmed by explicit cal- 
culation. Since the flux line displacement correlations 
< [u{x) — m(0)]2 > « £_q(x/Rl) in the Larkin regime and 
Ri [log{Ax / Rc) + B] in the logarithmic regime {A and B 
are constants [5]), 5 can be calculated exactly in those 
two regimes and the contributions of these two regimes 
are indeed negligible. 

We note that using these numerical values. 
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order of lO^T and the Bragg glass to vortex 
glass transition magnetic field Hm{T = 0) = 

{{7rcL)y{l6TT)^/V)){eo/UpfH!pH^lL, where Cl is 
the Lindeman criterion, is ^ 12T (for Hc2 = lOOT and 
Cl = 0.12) [5,6]. Similar values but Up/eo = 1/60 leads 
to Hm{T = 0) « 20T, and if 7 = 35 and Up/eo = 1/40 
one obtains Hm{T = 0) « 4r and Hcross « 18T (we 
have used these three different cases in Fig. 2). 

Focusing now on the RM regime c ontribution to 6, 
the displacement correlation function {[u{x) — w(0)]2) = 
iao/Trf{x/Rc)^'^/''\ where v = 1/6 [5,6,17]. The mag- 
netic field variance can then be calculated from Eq. (8): 



c2 ^1 

Orm -Z^^g 



,2 



^sm{XGy)dy, (11) 



where a = Rl/X, b = Rc/X, m = SG'^/{3K^). The 
behavior of the magnetic field variance (from numerical 
integration) is shown in Fig. 2. 

One notes that the magnetic field variance increases 
by increasing either magnetic field or the strength of the 
impurities. Since Rl is very small and Rc is very large, 
we can approximately replace the limit of integral from 
to 00, and do the integral analytically. The results 
of analytical integration show the same behavior for the 
magnetic field variance, in agreement with numerical in- 
tegration. 
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Fig. 2: Behavior of the variance S of the internal magnetic 
field as a function of the strength of the applied magnetic field 
for different strength of the disorder. 

We now briefly discuss our results in the context of 
/iSR measurements [9,11]. In Ref. [11], an increase of 
A°*^(5o) at T = 4 K between 25 A to 50 A per Tesla was 
reported. This corresponds to an effective narrowing of 
/iSR linewidth. Above, in Fig. 2, we found a broaden- 
ing of the internal magnetic field distribution by approx- 
imately 5—10 Gauss per Tesla of applied field Bq. To 
compare the magnitude of our field induced linewidth ef- 
fect to the experiment, we invoke the formula for S that 
applies to that of a perfect Abrikosov triangular vortex 
lattice in the London limit ^avl ~ O.OGl^a/ \^ , where 
$0 = 2 X IQ-^^ Tm^ [9]. We find that an increase of 
A''^ of between 25l to 50l reported in Fig. 3 of Ref. 
[11] corresponds to a decrease of the width of the inter- 
nal magnetic field distribution by about 4—8 Gauss per 
Tesla of applied field. In other words, the effect we report 
is of the same magnitude as that reported in /xSR exper- 
iments [11] but ascribed therein to the intrinsic nature 
of the field dependence of the microscopic A. The fact 
that we get a contribution of opposite trend (i.e. broad- 
ening) to that observed in experiments is irrelevant. Our 
results suggests that the extrinsic modification of the in- 
ternal field distribution as the system evolves from the 
Bragg glass to the vortex glass is of the same order as a 
contribution that may be of microscopic origin. Conse- 
quently, a quantitative description of the /xSR data that 
incorporates both intrinsic and extrinsic effects is desir- 
able in order to expose in a quantitative manner any 
microscopic effects such as those discussed in Refs. [2] 
and [3]. As mentioned above, the fiSK data are typically 
parametrized by including a phenomenological Gaussian 



broadening parameter [9-11]. Possibly, the phenomenol- 
ogy exposed in our calculation above may already be han- 
dled to some extend by the usage of a phenomenogical 
Gaussian broadening parameter. However, more theo- 
retical, numerical and experimental work is needed to 
assess whether or not this is the case. To conclude, we 
have shown that there can be significant magnetic field 
dependence for the magnetic field variance in the Bragg 
phase which is extrinsic and disorder driven and has no 
microscopic d— wave origin. We hope that our prelimi- 
nary study will motivate further work. 
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